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ABSTRACT 

We use ultradeep HST WFC3/IR imaging of the HUDF to investigate the rest-frame optical morphologies of 
a mass-selected sample of galaxies at z ~ 2. We find a large variety of galaxy morphologies, ranging from large, 
blue, disk-like galaxies to compact, red, early-type galaxies. We derive rest-frame u—g color profiles for these 
galaxies and show that most z ~ 2 galaxies in our sample have negative color gradients such that their cores 
are red. Although these color gradients may partly be caused by radial variations in dust content, they point 
to the existence of older stellar populations in the centers of z ~ 2 galaxies. This result is consistent with an 
"inside-out" scenario of galaxy growth. We find that the median color gradient is fairly constant with redshift: 
(A(w-g res t)/A(logr)) me di an = -0.46, -0.44 and -0.49 for z ~ 2, z ~ 1 and z = 0, respectively. Using structural 
parameters derived from surface brightness profiles we confirm that at z ~ 2 galaxy morphology correlates well 
with specific star formation rate. At the same mass, star forming galaxies have larger effective radii, bluer 
rest-frame u-g colors and lower Sersic indices than quiescent galaxies. These correlations are very similar to 
those at lower redshift, suggesting that the relations that give rise to the Hubble sequence at z = are already in 
place for massive galaxies at this early epoch. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: high- 
redshift 



1. INTRODUCTION 

Since the description of the Hubble sequence (e.g., B ubble 
Il926t Ide Vaucou leurs 1959; Sandage & Tammannl l 19811) we 
have learned that position along the Hubble sequence cor- 
relates wit h parameters like color, stellar age, and gas frac- 
tion (e.g., [Roberts & Ha vnesl 119941) . However, morpholo- 
gies by themselves provide very limited information, as they 
are scale free and do not include physical parameters such 
as surface brightnesses, sizes, luminosities and masses. Us - 
ing the Sloan Digital S ky Survey (SDSS, I York et al.ll2000b . 
Kauffman n et alj ( 120031) showed that the main parameter driv- 
ing galaxy properties is stellar mass. High mass galaxies are 
generally red, have old stellar populations and low specific 
star formation rates (SSFRs), while low mass galaxies are 
generally blue, have young stellar populations and high SS- 
FRs. 

A key question is what the structure was of the progeni- 
tors of low redshift galaxies. Hubble Space Telescope studies 
out to redshift z ~ 1 indicate that the morphologic al variatio n 
is comparable to that at low redshift (e.g., Bell et al. 2004). 
More recent studies of z > 2 galaxies using the HST NICMOS 
camera have yielded varying r esults, largely due to d ifferent 
selection criteria. For example, Pap ovich et al.l (120051) studied 
the rest-frame optical morphologies of a flux-limited sample 

1 Based on observations with the NASA/ESA Hubble Space Telescope, 
obtained at the Space Telescope Science Institute, which is operated by 
AURA, Inc., under NASA contract NAS 5-26555. 

2 Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Lei- 
den, The Netherlands 

3 Department of Astronomy, Yale University, New Haven, CT 06520- 
8101, USA 

4 Carnegie Observatories, Pasadena, CA 91101, USA 

5 UCO/Lick Observatory, University of California, Santa Cruz, CA 
95064, USA 

6 Center for Astrophysics and Space Astronomy, University of Col- 
orado, 389-UCB, Boulder, CO 80309, USA 



of galaxies at z ~ 2.3 and found that they are generally irregu- 
lar! jTbftetal] d2005l) . on the other hand, investigated the rest- 
frame optical and UV morphologies of distant red galaxies 
(DRGs) in the Hubble Ultra Deep Field (HUDF), and found 
both galaxies with irregular morphologies and galaxies with 
smooth morphologies. Additionally, they showed that the 
rest-frame optical morphologies of these galaxies are much 
more regular and centrally concentrated than the rest-frame 
UV morphologies. 

With the advent of the Wide Field Camera 3 (WFC3), with 
its vastly improved sensitivity and resolution compared to 
NICMOS, it has become possible to analyze the rest-frame 
optical structure of hi gh redshift galaxies w ith an unprece- 
dented level of detail. Came ron et al.l (|2010) have used data 
from the first year of observations of the HUDF and the Early 
Release Science Field to classify the rest-frame UV and op- 
tical morphologies of galaxies up to z ~ 3.5. These authors 
confirm results by e.g. iKriek et ail (120091) . who showed that 
massive galaxies at z ~ 2.3 can be separated into two dis- 
tinct classes: blue star-forming galaxies with irregular mor- 
phologies on the one hand, and red quiescent galaxies with 
smoother morphologies on the other. 

In this Letter, we extend the previous results using the full 
two-year ultradeep near-infrared (NIR) imaging of the HUDF 
taken with the HST WFC3. These data are the deepest ever 
obtained in the NIR and make it possible to analyze the mor- 
phologies, colors and structure of galaxies to z ~ 3 in the rest- 
frame optical. Using the incredible sensitivity and angular 
resolution of the WFC3 images we analyze the rest-frame op- 
tical surface brightness profiles of a mass-selected sample of 
galaxies at z ~ 2. We use these profiles to derive structural pa- 
rameters such as size and profile shape, and obtain rest-frame 
color profiles. We study the correlations between these pa- 
rameters as a function of redshift in order to investigate the 
Hubble sequence at different epochs in the history of the Uni- 
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FIG. 1 . — (a) Morphologies of z ~ 2 galaxies in the HUDF. The blue, green and red color channels are composed of PSF-matched HST/ACS (775, HST/WFC3 
F105 and HSTAVFC3 ffiso images, respectively. Mass increases upwards. The galaxies exhibit a very large range in morphologies: from very compact, nearly 
unresolved spheroids with red colors, to large, extended disk-like galaxies with blue colors. Most of the galaxies have a well-defined center, which is usually 
redder than the outer parts of the galaxy, (b) Four massive galaxies at z ~ 2. From left to right, the galaxies are shown as observed in the HST/ACS V506 band, 
the HST/WFC3 ffieo band, and the (ground-based) K s band. The difference between the (PSF matched) V6O6 band and //160 band images is large: whereas the 
galaxies exhibit very regular morphologies and have well-defined centers in the H[(,q band, they are very clumpy and irregular in the 1406 band, and in some 
cases are nearly undetected. These complex, wavelength-dependent morphologies suggest that z. ~ 2 galaxies may be composed of multiple components with 
very different stellar ages. 



verse. Throughout the Letter, we assume a ACDM cosmology 
with n m = 0.3, Q A = 0.7 and H = 70 km s" 1 Mpc" 1 . 

2. MORPHOLOGIES AT Z - 2 

We use the full two-year data taken with WFC3/IR of the 
HUDF, obtained in 2009-2010 as part of the HUDF09 HST 
Treasury program (GOl 1563). It consists of 24 orbits of F105 
imaging, 34 orbits of /125 imaging, and 53 orbits of Z/160 imag- 
ing. These images were reduced using an adapted pipeline 
dBouwens et al.ll2010l lOesch et al.ll20ia) . The FWHM of the 
point-spread function (PSF) is w 0.16 arcsec. The images 
are combined with very deep Advanced Camera for Surveys 
(ACS) images of the HUDF dBeckwith et alj 120061) to con- 
struct color images of the massive g alaxies in the field. 

We use the ^-selected catalog of [Wuvts et al.1 d2008l) to se- 
lect galaxies in the HUDF for which WFC3 imaging is avail- 
able. This catalog combines observations of the Chandra 
Deep Field South ranging from ground-based U band data 
to Spitzer 24/im data, and includes spectroscopic redshifts 
where avail able, as weU as pho tometric redshifts derived us- 
ing EAZY (Bramm er et al.l 120081) . Stellar masses were esti- 
mated from spectral energy distribution fits to the full pho- 
tometric data set (Forster Schreiber et al., 2011 in prepara- 
tion), assuming a Kr oupa initial mass function and the stellar 
population models of lBruzual & Ch ariot (2003). To study the 
morphological variation at z ~ 2, we select the 16 most mas- 



sive galaxies with 1.5 < z < 2.5. These galaxies have stellar 
masses between 1.2 x 10 10 M Q to 1.3 x 10 11 M Q . Color im- 
ages of these galaxies are shown in Figure QJ. A summary of 
their properties is given in Table [1] 

From Figure QJ it is apparent that galaxies at z ~ 2 show a 
large variation in morphology, size and color. One can dis- 
tinguish red, smooth, compact galaxies; blue galaxies with 
disk-like structures, some even with apparent spiral arms; and 
other star forming galaxies which appear more irregular. Mo st 
of the galaxies have a well-defined, red center. This is further 
illustrated in Figure [Tb, where galaxy morphology is shown 
as a function of wavelength. Four massi ve galaxies are show n 
in the observed HST/ACS V m band dBeckwith et alj|2006l) . 
HSTAVFC3 H\(,o band and K s band (based on ground-based 
imaging by Labbe et al., in preparation). The morphology 
of the sources differs strongly as a function of wavelength. 
The rest-frame UV morphology is always more clumpy and 
extended than the rest-frame optical morphology. In addition, 
the rest-frame optical images show well-defined centers for all 
our sources, whereas these are very often lacking in the rest- 
frame UV images. This confirms that the mass distribution 
of z ~ 2 galaxies is smoother and more centrally concentrated 
than would be concluded from rest-fram e UV imaging (see 
e.g. JLabbe et al.ll2003bllToft^tini2005h . 

3. COLOR GRADIENTS AT Z - 2 



Morphological evolution of galaxies to z ~ 2. 
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FIG. 2. — Top panels: surface density profiles of galaxies at z = 0, z ~ 1 and z ~ 2, color-coded according to specific star formation rate. The profiles are 
plotted up to the radius where uncertainties in the sky determination become significant. Surface density has been calculated from the average M/Lh (M/L r for 
the SDSS galaxies) of the galaxy multiplied by the observed H[(,q (r&s) band surface brightness profile, ignoring gradients in the mass-to-light ratio. It is clear 
that a large diversity of profiles exists, with some close to exponential profiles and others close to r 1 / 4 profiles. This diversity extends to z ~ 2. Middle panels: 
rest-frame u—g color profiles of the same galaxies. The profiles are normalized so that Au-j rfS = at r — r e . The color profiles are plotted from the PSF HWHM 
(~ 0.6 arcsec for the SDSS galaxies, ~ 0.08 arcsec for the HUDF galaxies) to the radius where the errors in the flux measurement reach 20%. The color profiles 
show an overall trend of redder colors at smaller radii, in all redshift intervals. Bottom panels: rest-frame u-g color gradients. Arrows indicate upper limits. 
Most z ~ 2 galaxies have negative color gradients. Color gradients do not seem to evolve very strongly between z ~ 2 and z = 0, for galaxies in the mass range 
considered here. 
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TABLE 1 
Galaxy properties. 



Source 


ID a 






n b 


'"stellar 


U-grest 


log 








(AB) 


(kpc) 




(Mo) 




(M© kpc" 2 ) 




HUDF09 


3203 


21.81 


3.49 


1.15 


10.65 


0.49 


8.77 


1.998 


HUDF09 


3239 


23.11 


0.32 


2.06 


10.51 


1.18 


10.72 


1.980 


HUDF09 


3242 


22.10 


0.44 


3.21 


10.76 


1.18 


10.67 


1.910 


HUDF09 


3254 


23.31 


1.96 


1.31 


10.25 


0.85 


8.87 


1.887* 


HUDF09 


3391 


23.28 


2.13 


0.44 


10.37 


-0.04 


8.92 


1.919* 


HUDF09 


3421 


23.79 


2.12 


0.90 


10.59 


0.19 


9.14 


2.457* 


HUDF09 


3463 


22.47 


1.80 


1.20 


10.09 


0.59 


8.78 


1.659* 


HUDF09 


3486 


22.43 


2.78 


0.65 


10.22 


0.88 


8.53 


1.628* 


HUDF09 


3595 


22.39 


3.67 


0.41 


10.98 


0.87 


9.06 


1.853* 


HUDF09 


3634 


22.94 


2.37 


1.28 


10.08 


0.58 


8.54 


1.994 


HUDF09 


3653 


23.85 


1.94 


0.95 


10.27 


1.24 


8.90 


1.776* 


HUDF09 


3721 


22.02 


3.00 


1.16 


10.54 


0.63 


8.79 


1.843 


HUDF09 


3757 


23.10 


1.25 


5.04 


10.28 


1.20 


9.29 


1.674* 


HUDF09 


3799 


24.59 


2.93 


1.50 


10.62 


0.87 


8.89 


2.492* 


HUDF09 


6161 


21.58 


6.70 


3.51 


11.08 


1.22 


8.63 


1.552 


HUDF09 


6225 


23.26 


2.58 


0.56 


10.33 


0.72 


8.71 


2.401* 


HUDF09 


6237 


23.51 


<0.10 


> 10.00 


10.72 


1.30 


14.40 


1.965* 



NOTE. — This is a sample of the full table, which is available in the online version of the article (see Supporting Information). 
a IDs correspond to NYU-VAGC IDs <Blanton et aI.I2005D for the SDSS galaxies, and FIREWORKS IDs tWuvts et ainoTjSl) for the HUDF galaxies 
b These parameters are derived from the SDSS g band imaging for the SDSS galaxies, and from the HST/W FC3 Hi^q band ima ging for the HUDF09 galaxies 
c Stellar masses and redshifts were obtained from the Guo et al. 1 2009) catalog for the SDSS galaxies, and the Wuyts et al. ( 2008) catalog for the HUDF galaxies. 
All masses include a correction factor to account for the difference between the catalog magnitude and our measured magnitude. 
d Save is defined as the average surface density within r e 
No spectroscopic redshifts are available for these galaxies; photometric redshifts are listed instead 
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In order to quantify the morphological properties of z ~ 2 
galaxies we measure their surface brightness profiles. In con- 
trast to conventional model-fitting techniques where a sim- 
ple model is used to approximate the intrinsic sur face bright- 
ness p rofile (e.g., using the GALFIT package of iPeng et al.l 
2002), we measure th e actual profile using the approach of 
Szom oru et al.l (12010b . The intrinsic profile is derived by fit- 
ting a Sersic model profile convolved with the PSF to the 
observed flux, and then adding the residuals from this fit to 
the unconvolved model profile. Effectively, the model pro- 
file is used to deconvolve the majority of the observed flux, 
after which this deconvolved profile is combined with the 
residuals to account f or deviations from the assumed model. 
iSzomoru et al.l (120101) have shown that it is thus possible to 
accurately measure the true intrinsic profiles out to large radii 
and very low surface brightness, even in cases where galaxies 
comprise a bright compact bulge and a faint extended disk. 

In addition to the z ~ 2 galaxies we measure the profiles of 
galaxies at z ~ 1 and z = 0. The z ~ 1 sample is taken from the 
same dataset as th e z ~ 2 sample, w hereas the z = galaxies 
are taken from the lGuo et all (12009b SDSS catalog of central 
galaxies. The galaxies are selected to lie in the same mass 
interval as the z ~ 2 sample: 1 .2 x 10 10 M Q < M stMar < 1 .3 x 
1O U M . The SDSS redshifts are required to be z = 0.03 ± 
0.015 and the z~ 1 galaxies have 0.5 < z < 1.5. This results 
in a sample of 27 galaxies at z ~ 1 and 84 galaxies at z = 0. 
Stellar masses and SSFRs for the S DSS sample are obtained 
from the MPA/JH U datareleaseQ (see B rinchmann et al.ll2004t 
ISalim et alJl2007l for details). A summary of the properties of 
the z ~ 1 and z = galaxies is given in Table Q] 

In Figure |2] we show surface density profiles, u-g rest color 
profiles, and u — g rest color gradients for galaxies in all three 
redshift bins. The surface density profiles are obtained by 
multiplying the surface brightness profiles in the H\eo band 
(g band for the SDSS galaxies) with the galaxies' average 
Msteiinr/ L ratios in that ban d (where M ste u ar is the total mass 
from the lWuyts et~aTll2008l catalog (MPA/JHU catalog for the 
SDSS galaxies) and L is derived from the surface brightness 
profiles). This approach ignores gradients in the M s , e u ar /L 
ratios and is therefore not exact; however, it allows for a 
more quantitative comparison between galaxies at different 
redshifts and with different colors. The color gradients are 
derived from fits to the u-g profiles between the half-width 
at half-maximum (HWHM) of the PSF (~ 0.6 arcsec for the 
SDSS galaxies, ~ 0.08 arcsec for the HUDF galaxies) and the 
radius where the errors in the flux measurement reach 20%. 

The majority of z ~ 2 galaxies in our sample have neg- 
ative color gradients, which do not vary strongly with red- 
shift: (A( M -g rest )/A(logr)) median = -0.47^, -0.33^ and 
-0.46^y 28 for z ~ 2, z ~ 1 and z = 0, respectively (where the 
errors give the 1-er interval around the median). Thus, galaxy 
color gradients seem to be remarkably constant with redshift, 
both for quiescent and star-forming galaxies. It should be 
noted that this is not a comparison of low-redshift galaxies to 
their high-redshift progenitors; the z ~ 2 galaxies are expected 
to evolve into more massive galaxies at low redshift, due to 
mergers and accretion. Since we use mass-limited samples, 
the z ~ 2 galaxies we consider will most li kely fall outside of 
our mass-limits at low redshift (see, e.g., Ivan Dokkum et all 
2010). 



http://www.mpa-garching.mpg.de/SDSS/ 



4. THE HUBBLE SEQUENCE FROM Z ~ 2 TO Z = 

Finally, we study the relations between structure, color and 
SSFR as a function of redshift in Figure [3] In the top row we 
show Sersic index against SSFR for the z = 0, z ~ 1 and z ~ 2 
samples. In the bottom row we show Sersic index against rest- 
frame u-g color for the same galaxies. The Sersic indices 
are derived from fits to the residual-corrected Hi6o-band (g- 
band for the z = galaxies) surface brightness profiles. There 
is a clear relation between these parameters at all redshifts: 
star-forming galaxies have "diskier" (lower n) profiles and 
bluer colors than quiescent galaxies. There is a large varia- 
tion in the SSFRs, colors and Sersic indices of z ~ 2 galax- 
ies, and the spread in these parameters is of roughly the same 
order of magnitude in all redshift bins. We note that the rela- 
tion shows systematic evolution: the median SSFR increases 
with increasing redshift and the median color and Sersic in- 
dex decrease, from log SSFR = -11.85^ yr" 1 at z = to 
log SSFR = -8.87!^" yr" 1 at z ~ 2, from u-g rest = 1.62^ 
to u-g rest = 0.87+g;g, and from n = 4.25^ to n = 1.20+§;^ 
(where the errors give the l-er interval around the median). 
These result s are qualita t ively consistent with the trends de- 
rived in e.g. JFranx et al.l (120081) . 

The relations between structure, color and SSFR are fur- 
ther illustrated in Figure |4] which shows rest-frame u-g 
color images of a selection of galaxies at z = 0, z ~ 1 and 
z ~ 2, as a function of SSFR. From each redshift bin we se- 
lect seven galaxies, evenly spaced in log SSFR. This figure 
very clearly illustrates the morphological variety present at 
Z ~ 2. Additionally, it confirms the relations shown in Fig- 
ure[3] star-forming galaxies are blue and extended, while qui- 
escent galaxies are red and relatively compact. Thus we find 
that the variation in galaxy structure at z ~ 2 is as large as at 
z = 0. Furthermore, the systematic relationships between dif- 
ferent structural parameters are very similar between z = and 
z ~ 2. This, in addition to the lack of evolution in the u—g rest 
color gradients shown in Figure |2] suggests that the underly- 
ing mechanisms that give rise to the Hubble sequence at z = 
may already be in place at z ~ 2. 

5. DISCUSSION 

We have shown that the morphologies of massive z ~ 2 
galaxies in the HUDF are complex and varied: from com- 
pact, apparently early-type galaxies to large star forming sys- 
tems superficially similar to nearby spirals. Many of these 
galaxies seem to be composed of multiple components with 
large differences in stellar age. We conclude that the vari- 
ety in morphologies which is observed at z = also exists in 
galaxies at z ~ 2. This is confirmed by an analysis of the 
surface density profiles, which reveals a large range in pro- 
file shapes. The profiles of the large star forming systems are 
close to exponential, whereas the profiles of the quiescent sys- 
tems are more concentrated. The correlations between mor- 
phology and SSFR are similar at all redshifts between z = 
and z ~ 2. 

This does not mean that morphologies and profiles are 
static: galaxies evolve quite strongly between z = 2 and z = 0. 
Quiescent galaxies are much more compact at high redshift, 
and cann ot evolve passively into low redshift quies c ent galax- 
ies (e. g.. iDaddi et all 120051; iTruiillo et all 120061: iToft et all 
120071: Ivan Dokkum et al.l 12008). A similar si ze evolution is 
requir ed for the star forming galaxies (e.g., [Williams et al.l 
2010). We note that due to the evolution in the mass func- 
tion of galaxies the z ~ 2 galaxies considered here will almost 
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FIG. 3. — Sersic index n plotted against SSFR (top row) and u-g rest color (bottom row) for galaxies with 1.2 X 10 10 Mq < Miliar < 1.3 X 1O H M0, at z = 0, 
Z ~ 1 and z ~ 2. Arrows indicate upper or lower limits. Profile shape and SSFR are anticorrelated in all three wavelength bins, whereas profile shape and color 
show a positive correlation. The relations between these morphological parameters are similar in all three redshift bins, although on average galaxies are bluer 
and have higher SSFRs and lower Sersic indices at high redshift. The similarity of the relations at z = and z ~ 2 suggests that the Hubble sequence was already 
in place for massive galaxies at z ~ 2. 
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FIG. 4. — Galaxy morphologies from z ~ 2 to z = 0, as a function of SSFR. The colors of the galaxies have been obtained by scaling their rest-frame u—g colors 
to the image's RGB space, as indicated by the color bar. The range in galaxy structure at z ~ 2 is very great, and is accompanied by a large range in galaxy colors. 
There is considerable evolution towards redder colors and lower SSFRs from z ~ 2 to z = 0. However, signs of a Hubble sequence (i.e., high SSFR galaxies are 
"diskier", more extended and bluer than low SSFR galaxies) appear to exist at z ~ 2. 
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certainly not evolve into the low-redshift galaxies in this pa- 
per. It would be very informative to examine the morpholo- 
gies of number density-limited galaxy samples over this red- 
shift range, in order to investigat e the evolution of the sam e 
galaxy population over time (e.g.. Ivan Dokkum et aTll20Tob . 

Additionally, we have derived u-g resl color gradients of 
our galaxy sample, and have shown that massive z ~ 2 galax- 
ies have negative color gradients that are comparable to those 
of low redshift galaxies in the same mass range. These color 
gradients may partly be caused by radial variations in dust 
extinction, but it is likely that radial changes in st ellar popu- 
lations play a large role (e.g., Abr aham et al.l ll999). This sup- 
ports a galaxy growth scenario where small galaxies formed 
at high redshift grow by accreting material onto their outer re- 
gions. However, more detailed studies are needed in order to 
disentangle dust and stellar age effects on the color gradients 
of high redshift galaxies. 

The results seem to contrast with ea rlier analyses of 
galaxies at z ~ 2 in the HDFN (e.g., iDickinsonl 120001 
Papov ich et al.| [2005). These authors concluded that all galax- 
ies at this redshift are irregular and compact, with little dif- 
ference between the rest-frame UV and optical . Field-to- 
field va riations may play a role; as shown in lLabbe et"aD 
(2003a) the HDFN contains very few massive high-redshift 
galaxies. Indeed, the galaxies studied by these authors have 
lower masses (M ste u ar < 3 x 10 10 M Q ) than the galaxies we 
consider, which could result in the difference with the study 
presented here. Additionally, We note that recent kinemati- 
cal studies of massive z ~ 2 galaxies als o indicate that many 
contain gas with ordered motion (e.g.. iGenzel et al.l 120081 
iForster Schreiber et al .1200 9, Tacco ni et al.l2010l) . The results 
presented in this Letter are fully consistent with those. 



Our results raise the question at what redshift the first "or- 
dered" galaxies appeared, with structures similar to the Hub- 
ble sequence. There are indications that at redshifts be- 
yond 3 such galaxies may be much harder to find. Typ- 
ical high redshift L yman break galaxies are very clumpy 
and irregular (e.g., lLowenthal et al.l 119971) . and differ sig- 
nificantly in appearance from regular spiral galaxies. Fur- 
thermore, the population of massive galaxies that are faint 
in the UV may be very small at redshifts beyond 3 (e.g., 
iBrammer & van Dok kum 2007]), at l east in the mass regim e 
that is considered in this paper (e.g., March esini et al.l 2010). 
With current observational capabilities we are severely lim- 
ited in studying the rest-frame optical properties of galaxies 
at redshifts beyond z ~ 3.5, due to rest-frame optical emission 
moving redward of the observers' K band. With improved ca- 
pabilities it may become possible to study the red massive 
galaxy population at these redshifts. Several candidates have 
been found at z > 5 with significant B aimer discontinuities 
(lEvles et alJl2005t iMobasher et al]l2005t iWiklind et al.ll2008t 
lRichardetal.ll201lb . and many more are speculated to ex- 
ist. These galaxies could very well be the centers of multi- 
component galaxies at redshifts between z = 3 and z = 5. 
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